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ABSTRACT 

We investigate detailed kinematical properties of simulated collisionless 
merger remnants of disk galaxies with mass ratios of 1:1 and 3:1. The sim- 
ulations were performed by direct summation using the new special hardware 
device GRAPE-5. In agreement with observations, the shape of the hne-of-sight 
velocity distribution (LOSVD) is Gaussian with small deviations. For most cases 
we find that the retrograde wings of the LOSVD are steeper than the prograde 
ones. This is in contradiction with observations which show broad retrograde 
and steep prograde wings. This serious problem in the collisionless formation 
scenario of massive elliptical galaxies can be solved if all rotating ellipticals, even 
boxy ones, contain an additional stellar disk component with oc 15% of the total 
stellar mass and a scale length of order the effective radius of the spheroid. We 
propose that the progenitor galaxies of massive ellipticals must have contained a 
significant amount of gas that did not condense into stars during the merger pro- 
cess itself but formed an extended gaseous disk before the star formation epoch. 
The heating source that prevented the gas from forming stars early and the origin 
of the large specific angular momentum required for the gas component to form 
an extended disk are still unsolved problems. 

Subject headings: galaxies: elliptical - galaxies: interaction- galaxies: structure 
- galaxies: evolution - methods: numerical 

1. Introduction 

The merger hypothesis (Toomre & Toomre 1972) assumes that elliptical galaxies form 
by major mergers of spiral galaxies. This is supported by observations of nearby merging 
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systems that seem to evolve into systems with properties comparable to elliptical galaxies 
(Hibbard & Yun 1999). The gravitational interaction during the merger process also 
triggers massive star formation and, probably, the formation of massive central black holes 
and AGNs (Sanders & Mirabel 1996; Rigopoulou et al. 1999). On the theoretical side, 
numerous numerical simulations have shown that mergers lead to spheroidal systems with 
global properties comparable to elliptical galaxies (Barnes & Hernquist 1992; Barnes 1998; 
Bekki & Shioya 1997; Naab, Burkert & Hernquist 1999). During the merging epoch the 
galaxies are far from dynamical equilibrium. Gas, if added to the simulations, condenses 
into stars by massive bursts of star formation during the merging process or effectively 
looses angular momentum and falls into the central parts of the merger remnant, probably 
feeding a central black hole (Mihos & Hernquist 1996; Springel 2000; Bekki 2000). When 
the merger remnants have settled into dynamical equilibrium, phase mixing and violent 
relaxation have erased most information about the initial conditions. However, as violent 
relaxation is incomplete, the stellar phase space distribution of elliptical galaxies should 
still show fine structure which affects the isophotal shapes and velocity distributions. This 
fine structure provides fundamental insight into the formation epoch of elliptical galaxies. 

Numerical simulations of coUisionless mergers of disk galaxies have confirmed this 
conclusion. The LOSVDs of merger remnants and the orbit structure of their stars contain 
information about the initial disk orientations (Heyl, Hernquist & Spergel 1996). The 
isophotal shapes and the dynamics of the merger remnants are primarily determined by 
their initial mass ratios (Naab, Burkert & Hernquist 1999; Bendo & Barnes 2000). 

In this letter we show that the LOSVDs of most eUiptical galaxies contain evidence 
for a second disk component with roughly 10% to 20% the luminosity of the spheroid 
and a large scale radius which if is order the effective radius of the spheroid. This disk 
might have formed by late gas infall from gas-rich progenitors, followed by star formation. 
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Observationally, photometric and kinematic investigations of elliptical galaxies have already 
indicated the existence of large stellar disks of the same order of magnitude as predicted 
here, especially in disky, rotating ellipticals (Rix & White 1990; Scorza et al. 1998). 
We demonstrate that this component is required in almost all ellipticals, independent of 
whether they are disky or boxy in order to bring the observations in agreement with the 
merger scenario of the formation of early-type galaxies. 

2. The merger models 

We use the method described by Hernquist (1993) to construct disk galaxies in 
dynamical equilibrium. The system of units is: gravitational constant G=l, exponential 
scale length of the larger disk in the merger h=l and mass of the larger disk — 1. The 
disks are exponential with an additional spherical, non- rotating bulge with mass M5 = 1/3, 
a Hernquist density profile (Hernquist 1990) and a scale length — 0.2h. The system lives 
in a spherical pseudo-isothermal halo with a mass — 5.8, cut-off radius Tc = lOh and 
core radius 7 = 1/1. 

The N-body simulations were performed using a direct summation code with the new 
special purpose hardware GRAPE-5 (Kawai et al. 2000). The 1:1 merger was calculated 
adopting in total 400000 particles with each galaxy consisting of 60000 disk particles, 20000 
disk particles and 120000 halo particles. For the 3:1 merger the parameters of the more 
massive galaxy were as described above. The low- mass galaxy contained 1/3 the mass and 
number of particles of the larger galaxy, with a disk scale length oi h = 1^1/3, as expected 
from the TuUy- Fisher relation. 

For both mass ratios, the galaxies approached each other on nearly parabolic orbits 
with an initial separation of 30 length units and a pericenter distance of 2 length units. In 
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this letter we focus on one corotating and one counterrotating geometry. We have found 
that these two impact parameters are representative for a much larger set of simulations 
that have been performed and that will be discussed in a subsequent paper. The inclinations 
of the corotating disk models (called model A for the equal mass merger and model C for 
the 3:1 merger) relative to the orbital plane are ti — 30° and t2 — —30° with arguments 
of pericenter of uji = 30° and 002 = —30°, respectively. The counterrotating disks (called 
model B for 1:1 and model D for 3:1) have h = 30° and ^2 = 150°, cui = 30° and U2 = -30°, 
respectively. In all simulations the merger remnants were allowed to settle into dynamical 
equilibrium for approximately 10 dynamical timescales after the merger was complete. 
Then their equilibrium state was analysed. 

3. LOSVD analysis of the merger remnants 

To measure the line-of-sight velocity distribution of a merger remnant we shifted 
the densest region of every two-dimensional projection to the origin. Then we placed 
a slit with a width of 0.2 unit lengths along the apparent long axis of each projected 
remnant. Thereafter we binned all particles falling within each grid cell in velocity along 
the line-of-sight. The grid spacing was chosen to be 0.15 times the projected half-mass 
radius. The width of the velocity bins was set to a value of 0.2 for hne-of-sight velocities v\os 
in the range —4 > Vi^s > 4. This results in 80 velocity bins over the whole velocity interval. 
Using the binned velocity data we constructed line-of-sight velocity profiles (LOSVD) for 
each bin along the grid. Subsequently we parametrized deviations from the Gaussian shape 
of the velocity profile using Gauss-Hermite basis functions (van der Marel & Pranx 1993; 
Gerhard 1993; Bendo & Barnes 2000). The kinematic parameters of each profile (cTfit, ffit, 
H3, H4) were then determined by least squares fitting. The large number of simulated 
stellar particles (> 100000) guaranteed that at least 2000 particles fall within each slit 
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inside one effective radius. 

4. Compcirison with observations 

The straight hnes in Figure 1 indicate the observed local correlation between and 
v/a for a sample of elliptical galaxies in low density environments, published by Bender, 
Saglia & Gerhard (1994). In all cases and v/a have opposite signs. The data indicate 
that all eUiptical galaxies have LOSVDs with steep prograde wings and broad retrograde 
wings. Mehlert et al. (2000) have shown that this result also holds for cluster ellipticals in 
Coma. 

The local correlation between H3 and v/a for the simulated corotating (A,C) and 
counterrotating (B,D) 1:1 and 3:1 mergers, respectively, is shown by the dots in Figure 
1. Every remnant is analysed as seen from 50 random viewing angles. For the corotating 
mergers the correlation between H3 and v/a is almost opposite to the observed one. This 
is also reflected in a positive effective value for (left plot in Figure 2) for almost all 
projection angles, where < > is defined as the mean value between 0.25 and 0.75 
effective radii (Bender, Saglia & Gerhard 1994). Here, a positive < > corresponds to an 
anticorrelation between H3 and v/a, such that i/3 > for f/cr < and vice versa. We find 
that < H^> versus v/a^ does not follow the observed correlation, indicated by the solid and 
dashed lines. The profiles have broad prograde wings and narrow retrograde wings. The 
counterrotating 3:1 merger (D) in Fig. 1 shows a large spread around zero but also does not 
agree with observations. The only exception is the equal mass merger of counter-rotating 
disks (B) which leads to a very anisotropic elliptical with no signature of rotation. 

Note, that the four simulations which we focus on in this paper have been chosen 
because they are representative for a much larger set of simulations with different orbital 
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geometries. These four models have been recalculated with high resolution for better 
statistics. 

We conclude that coUisionless mergers of disk galaxies in general fail to explain the 
detailed kinematics of all observed elliptical galaxies (also massive, boxy ones) that exhibit 
a significant amount of rotation {y/a > 0.2) inside one r^s- Further evidence for a possible 
failure of the coUisionless merger picture comes from Cretton et al. (2001) who showed that 
the kinematical properties of very faint and fast rotating, disky elliptical galaxies can not 
be explained by coUisionless mergers of disk galaxies. 

5. Theoretical evidence for disks in elliptical galcLxies? 

One plausible explanation for line-of-sight velocity distributions with negative is the 
superposition of a spheroidal body with a disk-like component (Bender, Saglia & Gerhard 
1994). With a simple experiment one can test if the wrong correlations found in our 
remnants result from a lack of such a component. We artificially added a thin (scale height 
~ O.OSreif), cold (velocity dispersion perpendicular to the disk — 0), stellar disk with an 
exponential surface density profile. The disk was placed in the plane defined by the long 
and intermediate axes of the main stellar body, rotating in the same direction as the main 
stellar body. The particles of the disk were assumed to move around the center of the galaxy 
in centrifugal equilibrium with the gravitational potential arising from the total enclosed 
mass. No additional random motion was added. Under these simple assumptions only two 
parameters remain free: the total mass of the disk Md in units of the total luminous mass 
of the remnant and its scale length in units of the projected half-light radius res- The 
results for the corotating 3:1 remnant (C) are summarized in the Figures 3 and 4. This case 
is representative for all our unequal mass merger remnants. In general, we find that disks 
with small masses or radii do not change the LOSVDs of the stellar component. Disks with 
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masses and radii in the region indicated by the black dots in Figure 4, on the other hand, 
lead to a significant change in the resulting line profile. The prograde wings become steeper 
than the retrograde ones in very good agreement with observations. This effect is shown in 
Figure 3 for a disk with 15% of the spheroid mass and = 1.25res- The infiuence on the 
< H3 > versus v/ao correlation and its dependence on projection effects is shown on the 
right hand side of Figure 2. Now the values fall in the observed regime. 

If the additional disk becomes too massive, the absolute values for are larger than 
observed and the surface brightness profiles change from de Vaucouleurs to exponential 
profiles which is again not in agreement with observations. We therefore conclude that the 
existence of an additional stellar disk component with 10% to 20% the luminosity of the 
spheroid and a scale length of the order of Tefr can explain the observed correlation between 
and v/a in elliptical galaxies. 

6. Discussion and Conclusions 

Our simulations indicate that models of pure coUisionless mergers in general fail to 
reproduce LOSVDs of observed eUiptical galaxies. All simulated profiles locally have the 
wrong sign of with respect to (v/a), compared with the observations. A similar result 
has been reported by Bendo & Barnes (2000). They, however, did not investigate projection 
effects. A detailed investigation of local properties of the LOSVDs and projection effects 
shows that the presence of an additional stellar disk component added to the stellar body 
after the merger is complete could solve the problem. 

One possible origin for such a disk is gas that must have been present during the 
formation epoch of massive ellipticals and that formed a disk after the major merger event 
was completed. In order to fit the observed profiles, the amount of gas that settled into 
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the equatorial plane must have been significant. It is puzzling why the gas did not turn 
into stars prior to its infall into the equatorial plane and how the gas could keep its angular 
momentum required for the disk to have a scale length equal to the effective radius of the 
spheroid. 

Barnes & Hernquist (1996) do find gaseous disk-like components in their equal mass 
merger models. Tidal torques lead however to efficient angular momentum loss in the 
gaseous component resulting in gas infall to the center. Only 20% of the initial gas mass 
(which was 10% of the initial disk mass) and by this less than 2% of the stellar mass settled 
into an extended disk- like component. In this case the disk would not be massive enough 
to change the LOSVD of the system significantly. Equal mass mergers of very gas rich 
disks (more than 50% gas) could, however, lead to the formation of more massive disk-like 
components in the end. 

First simulations by Naab, & Burkert (2001) have shown that for unequal mass 
mergers of gas rich galaxies which lead to disky, fast rotating ellipticals due to the stronger 
centrifugal support, a large fraction of the gas settles into extended disks after the merger 
is complete, if star formation is suppressed. This seems to be an attractive scenario to 
explain the existence of large disks in disky elliptical galaxies. The origin of large stellar 
disks in massive boxy ellipticals with a small but significant amount of rotation that form 
presumably from equal-mass mergers is, however, still unclear. 



We thank Ralf Bender, Hans- Walter Rix and Nicola Cretton for helpful discussions. 
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Fig. 1. — Local correlation between and v/a for the 1:1 and 3:1 merger simulations. The 
models (A,C) and (B,D) have co -and counterrotating geometries, respectively. The dots 
represent the data for each model as seen from 50 random viewing angles. The error bar in 
the upper panel shows the bootstrap error at 0.75 Vcs- The observed correlation from Bender, 
Saglia & Gerhard (1994) is tentatively indicated by the straight lines. For the complete data 
set see their Figure 15. 
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Fig. 2. — Left: Effective versus v/ao (rotational velocity at one res over central velocity 
dispersion) for 50 random projection of model (C). Right: Same for model (C) with an 
additional disk as described in the text. The straight line is the mean correlation observed 
by Bender, Sagha & Gerhard (1994). The parallel long-dashed lines give an indication of 
the observed spread. The short dashed line is < i^a > = 0. 



-14- 




Fig. 3. — Local correlation between H^, and v/a for model (C) with an additional thin stellar 
exponential disk with a scale length of 1.25refr. A mass of 15% of the spheroid is assumed. 
The error is derived in the same way as in Figure 1. 
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Fig. 4. — Mass versus scale radius for a disk added to the merger remnant (C). 
Combinations that are able to reproduce (big dot), almost reproduce (small dot) or fail to 
reproduce (minus sign) the observed correlations are shown. 



